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Abstract 

The expOU stochastic volatility model is capable of reproducing fairly well most important 
statistical properties of financial markets daily data. Among them, the presence of multiple time 
scales in the volatility autocorrelation is perhaps the most relevant which makes appear fat tails in 
the return distributions. This paper wants to go further on with the expOU model we have studied 
in Ref. by exploring an aspect of practical interest. Having as a benchmark the parameters 
estimated from the Dow Jones daily data, we want to compute the price for the European option. 
This is actually done by Monte Carlo, running a large number of simulations. Our main interest 
is to "see" the effects of a long-range market memory from our expOU model in its subsequent 
European call option. We pay attention to the effects of the existence of a broad range of time 
scales in the volatility. We find that a richer set of time scales brings to a higher price of the option. 
This appears in clear contrast to the presence of memory in the price itself which makes the price 
of the option cheaper. 
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I. INTRODUCTION 



The model, suggested by Bachelier in 1900 as an ordinary random walk and redefined 
in its final version by Osborne in 1959 j^, presupposes a constant "volatility" a, that is 
to say, a constant diffusion coefficient D = a 2 . However, and especially after the 1987 
crash, there seems to be empirical evidence, embodied in the so-called "stylized facts" , that 
the assumption of constant volatility does not properly account for important features of 
markets 0,0,0]. It is not a deterministic function of time either (as might be inferred by 
the evidence of non stationarity in financial time series) but a random variable. In its more 
general form one therefore may assume that the volatility a is a given function of a random 
process Y(t), i.e., a{t) = (p(Y{t)). 

At late eighties different stochastic volatity (SV) models [Brownian motion with random 
diffusion coefficient] were presented for giving a better price to the options somewhat ignoring 
their ability to reproduce the real price time series [6|. More recently SV models have 
been suggested by some physicists as good candidates to account for the so-called stylized 
facts of speculative markets ^,0, 10, 11, 12, In the context of mathematical finance, 
we mainly have three models: the Ornstein-Uhlenbeck (OU) [llj], the Heston and the 
exponential Ornstein-Uhlenbeck model [jfl. We have recently tried to decide which model 
works better 0, Q|. Very recently we have studied the exponential Ornstein-Uhlenbeck 
stochastic volatility model Q| and observed that the model shows a multiscale behavior in 
the volatility autocorrelation [jj. It also exhibits a leverage correlation and a probability 
profile for the stationary volatility which are consistent with market observations. All these 
features seem to make the model more complete than other stochastic volatility models also 
based on a two-dimensional diffusion. It is worth to mention that there has been some more 
sophisiticated models based on a three-dimensional diffusion process that reproduce the 
intrincate set of memories in market dynamics Q,^^. Indeed, coming from multifractality 
framework, there are recent nice papers with promising features to provide even a better 
description that the one by the common stochastic volatility models |18| . 
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II. THE VOLATILITY MODEL 



Let us briefly summarize the main definitions and the properties of the exponential 
Ornstein-Uhlenbeck stochastic volatility model [lj. The model consists in a two-dimensional 
diffusion process given by the following pair of Ito stochastic differential equations (SDEs): 

X(t) = me Y ^i(t) (1) 
Y(t) = -aY(t) + k£ 2 (t), (2) 

where dots represent time derivative. The variable X(t) is the undrifted log-price or zero- 
mean return defined as X = S(t)/S(t) — (S(t)/S(t)\ where S(t) is a the asset price. The 
parameters a, m, and k appearing in Eq. (J2J are positive and nonrandom quantities. The 
two noise sources of the process are correlated Wiener processes, i.e., are zero- mean 
Gaussian white noise processes with cross correlations given by 

m)ut))=PiAt-t% (3) 

where pa — 1, pij — p (i ^ j, — 1 < p < 1). In terms of the proces Y(t) the volatility is given 
by 

a(t) = me Yit) . (4) 



It is worth to mention that multifracta 
describes the logarithm of the volatility 



s models also considers a random variable which 



3, lis i . 



Among the most important results of the expOU, we must stress the following three Q . 
First of all the stationary volatility pdf which is a log-normal distribution being quite con- 
sistent with empirical data p| 

Pst{<r) = ^=exp{-ln 2 (a/m)/2/3}, (5) 

where (3 = k 2 /2a. Notice that the stationary distribution broadens the tails as we increase 
the value of (3. Secondly, we have the leverage effect 

£(r) = (dX fd7m 2T)2) = {2pk/m) exp { _ar + me ~ aT ~ 3/4) ) ^ (r) ' (6) 

where H(r) is the Heaviside step function. One can show that 

C(t) ~ -pkeP /2 e- k ' 2r (or < 1), (7) 
m 



FIG. 1: Volatility autocorrelation (JSJ) for different values of (3. The characteristic time scale for 
the long time corresponds to 1/q while short time scale in leverage is related to 1/k 2 . For a large 
value of /3 we observe a richer multiple time scale behaviour. 

while leverage is negligible for large times. These approximations hold only if we take /3 > 1. 
Thirdly, and perhaps the most important, the volatility autocorrelation 

(dX(t) 2 dX(t + r) 2 )-(dX(t) 2 } 2 = exp[4/3e— ] - 1 
lTj (dX{ty) - (dX{t) 2 ) 2 3e^-l ' { ) 

which expanded in the right way allow us to observe a cascade of exponentials 

Co^H^g^e— . (9) 

This expression indicates that there are infinite time scales between the limiting behaviours 

Corr(r) « [or » 1), Corr(r) = [e 4 ^ 2fcV - l] + (a 2 r 2 ) 

(10) 

As one can observe the characteristic time scale for the long time might be associated to 
1/a while short time scale in leverage is related to 1/k 2 (see Fig. Q). The distance between 
the smallest time scale and the largest is given by (3 = k 2 /2a > 1. The bigger j3, the larger 
is the distance and the richer is the cascade of multiple time scales of the process. Even 
more, as we enlarge the distance between smaller time scale and larger time scale, we also 
get a broader lognormal distribution (cf. Eq. (J3J)) for the volatility and a fatter distribution 
fo the returns. 
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III. OPTION PRICING 



Having in mind latter expressions and main messages concerning (3, we can start looking 
at the inferrences of the volatility long range memory in the option price. An European 
option is a financial instrument giving to its owner the right but not the obligation to buy 
(European call) or to sell (European put) a share at a fixed future date, the maturity time 
T, and at a certain price called exercise or striking price K. In fact, this is the most simple 
of a large variety of derivatives contracts. In a certain sense, options are a security for the 
investor thus avoiding the unpredictable consequences of operating with risky speculative 
stocks. 

The payoff of the European call contract at maturity date T is 



To compute the price of the call we can use the average (on an equivalent martingale measure) 
over a large set of simulated paths which can be written mathematically as follows 



To do so we assume the process defined in the pair of Eqs. (JU)-© but with a drift equal to 
the risk free interest ratio r and with a market price of risk for the volatility component set 
to be 0. That is: keeping the current measure for the dynamics of the volatility. There are 
many subtleties behind the evaluation of the volatility market price of risk and we address 
the reader to the book by Fouque et al. to get a deeper knowledge on this issue. 

We thus simulate the random path assuming one day timesteps. Figure El shows the 
different paths observed for the SV model and the log-Brownian one using m as a volatility 
o diffusion coefficient (cf. Eq.(0J)). Repeating the paths 1,000,000 times over the same time 
window, we can get an approximate idea of the correction to the classic Black-Scholes price 
formula. We recall that Black-Scholes assumes a log-Brownian motion for the underlying 
whose price is well-known and has the analytical solution: 



S(T) - K if S(T) > K, 



and otherwise. 



C T (S) = E* e- rT (S(T) - K) + S(t = 0) = S . 



(11) 



C BS (S,t) = SNidJ 



Ke- r ^N{d 2 ) (0 < t < T), 



(12) 



where N(z) is the probability integral 




(13) 
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FIG. 2: Random path of the log-Brownian motion compared to a couple of expOU SV model 
paths with different /3's. Fluctuations increases wildly with (5. Price time series take returns with 
a non-zero drift r = 1.845 day -1 (5% annual risk free interest rate ratio). Rest of the parameters 
are k = 0.18 day _1//2 , m = 0.00141 day~ 1//2 and for the log-Brownian motion (constant volatility) 
we take a = m. All series take an initial stock price S = 0.999-fT. 



and its arguments are 



di 



ln(S/K) + (r + a 2 /2)(T - 1) 



d 2 = d x - ay/T - t. 



(14) 



In contrast with some contributions from mathematical finance, we are not inserting the 
parameters values blindly nor providing a large collection of parameters where it is quite 
hard to intuit the meaning and the effects of each parameters. We take the parameters set 
in Ref. for the Dow Jones daily index as a benchmark. The parameters derived gives an 
opposite approach to the one already performed by Fouque etal. [t| for the expOU. They are 
focused in analytical asymptotic results with the cases where (3 < 1 but the problem is that 
with this restriction one does not have the desired cascade of time scales in the volatility 
aut o correlat ion . 

We recall that m is the volatility level assuming no auxilary diffusion process for the 
volatility, p gives the asymetry between positive and negative return change and the appro- 
priate sign for the leverage effect. And finally, short range memory time scale is measured 
by k 2 which appears to be of the order of few weeks (10-20 days). We will now focus on the 
effect of the largest time scale memory a -1 which is considered to be of the order of years. 
We have shown in Ref. [lj that a good approximate solution for reproduce the memories of 
the market is performed taking k 2 = 1.4 x 10~ 2 day -1 , p = 0.4 and m = 1.5 x 10~ 5 day -1 / 2 . 
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FIG. 3: The European call option in terms of the moneyness S/K. We show the Black-Scholes 
formula and the expOU SV model for several values of the parameters. The computation with the 
expOU is performed by running 1,000,000 paths and average over their final premium. Results 
become very noisy for (3 around 10. 

The comparison between shortest and largest time scale is provided with (3 and at least for 
the Dow Jones daily data this is around 3.8. 

In FigureHHwe plot the call price for several values of (3 averaging over the initial volatility 
with the lognormal distribution (|5|). We take maturity date at 20 days and we represent the 
option assuming that the risk free interest ratio to 5% per year. Even for such a small time 
horizon (much smaller than the volatility memory), we get important price differences. And 
in any case the longer the memory is the more expensive is the option. This can also be 
quantified by the relative price between new option price and the Black-Scholes one, that is: 

C{S,t) 1 



Relative difference 



(15) 



C B s(S,t) 

The results are plotted in Fig. One observe that the price difference becomes more 
important with a monotonic increase in terms of (3. These differences may become relatively 
very important for small moneyness S/ K. And the decay of the relative difference for larger 
distances with respect to the striking price K is becoming slowler with a higher value of /3. 
We have tested the results with different maturity dates and with different values for k, p, 
and m with similar conclusions. 

In a previous paper 19] , we also have been studying the effects of the memory in the option 
price. In that case, however, we had memory in the price itself breaking the efficient market 
hypothesos. We had observed that the call became cheaper with this kind of memory even 
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FIG. 4: The relative difference between the expOU call option and the European call option in 
terms of the moneyness S/K. We show the difference for several values of the parameters. The 
computation with the expOU is performed by running 1,000,000 paths and average over their final 
premium. 

if this of only one day. The presence of memory in the volatility has opposite consequences. 
This paper has aimed to insist in the fact that the memory persistence in volatility affects 
the price making this to be higher. 
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